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ABSTRACT 

Radio interferometry is a powerful technique for astronomi- 
cal imaging. The theory of compressed sensing (CS) has been 
appHed recently to the ill-posed inverse problem of recover- 
ing images from the measurements taken by radio interfero- 
metric telescopes. We review novel CS radio interferometric 
imaging techniques, both at the level of acquisition and re- 
construction, and discuss their superior performance relative 
to traditional approaches. In order to remain as close to the 
theory of CS as possible, these techniques necessarily con- 
sider idealised interferometric configurations. To realise the 
enhancement in quality provided by these novel techniques 
on real radio interferometric observations, their extension to 
realistic interferometric configurations is now of considerable 
importance. We also chart the future direction of research re- 
quired to achieve this goal. 

Index Terms — Compressed sensing, sparsity, radio 
interferometry, interferometric imaging 

1. INTRODUCTION 

Aperture synthesis interferometry is a powerful technique in 
radio astronomy, providing astronomical images at angular 
resolution and sensitivity that are otherwise inaccessible (see 
e.g. 1 1 1). Arrays of telescopes track a source as it traverses the 
sky to synthesise the aperture of a telescope of much larger 
size than the individual telescopes comprising the array. Ra- 
dio interferometric arrays directly probe the Fourier plane as- 
sociated with the image tangent to the celestial sphere at the 
source position. As the source traverses the sky, interfero- 
metric measurements, termed visibilities, are generated on an 
elliptical track of positions in the Fourier plane. Visibility 
coverage in the Fourier plane is therefore incomplete. Con- 
sequently, recovering interferometric images from measured 
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visibilities requires solving an ill-posed deconvolution prob- 
lem. 

The recovery of interferometric images has traditionally 
been tackled by the so-called CLEAN algorithm |2|, which 
is a local matching pursuit approach. Recently, however, the 
theory of compressed sensing (CS) U IH has been applied 
to interferometric imaging with considerable success, both at 
level of acquisition and reconstruction [ 5^,^,71. Furthermore, 
CS techniques have been developed to successfully extract 
astronomical signals of interest from interferometric observa- 
tions corrupted by background contributions 1 8 1 . First steps 
in these directions necessarily involved the study of idealised 
interferometric configurations to remain as close to the theory 
of CS as possible. With a new generation of interferometric 
telescopes under design and construction, radio interferom- 
etry will continue to play a fundamental role in astronomy. 
The extension of CS techniques to more realistic interfero- 
metric configurations is therefore now of paramount impor- 
tance, in order to realise the enhancement in quality that these 
techniques provide on real observations made by radio inter- 
ferometric telescopes. 

In this article we review the current state-of-the-art of CS 
imaging techniques for radio interferometry. We then chart 
the future direction that this research must take to be appli- 
cable to real radio interferometric telescopes. In Sec. [2] we 
review radio interferometry concisely and pose the interfer- 
ometric imaging inverse problem. We discuss the CS tech- 
niques applied to this problem in Sec. |3j before charting the 
future direction of research in Sec.|4l 

2. INTERFEROMETRIC INVERSE PROBLEM 

The field-of-view (FoV) observed by an interferometer is lim- 
ited by the primary beam of the telescope A. We denote the 
sky intensity to be imaged by x and consider a coordinate 
system centred on the pointing direction of the telescope Sq. 
where a point s on the celestial sphere is defined by its co- 
ordinates {l^m^n), noting n = n{l) = (1 — H^Hi)"^^^. for 
/ = (/,m) and ||/||2 = + m^. This configuration is illus- 
trated in Fig.[T] The complex visibility measured by the cor- 
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Fig. 1. Coordinate system for interferometric imaging (credit: 
|6)). 

relation of the incoming electric fields of the observed source 
for each telescope pair is given by 111 

y{uM=l A(0x(0C(||/||2)e-^2-^-^^, (1) 

where the measured visibility depends only on the rela- 
tive positions between telescope pairs, denoted in units of 
wavelengths by the baseline vector 6a = (u^v^w), with 
u = (li, v), and the so-called -modulation is given by 

Visibilities are acquired for non-zero w if the baseline of the 
telescope pair has a component in the pointing direction of 
the telescope. 

Typically small FoV assumptions are made, in which case 
the approximations n{l) 1 and C(||/||2) — 1 reduce ([T]) to 
a standard Fourier transform. However, if these assumptions 
are relaxed, n~^{l) may instead be absorbed into the beam 
but a -modulation remains. The i(; -modulation gives rise to 
the spread spectrum phenomenon, which is a serendipitous 
property of the acquisition, enhancing the performance of CS 
techniques considerably. 

To recover the source image from incomplete visibility 
measurements, we pose the inverse problem ([T]) in a discrete 
setting. The following formulation is restricted to constant w, 
thus restricting visibilities to a Fourier plane We measure M 
discrete visibilities denoted by the vector y, falling on a grid 
of frequencies of size A^ = 7V^/^xA^^/^, which are related to 
the N source image pixels, denoted by the vector x, through 
the linear system 

y = ^x^n, (2) 

with 

<!> = MFCA, 

^The restriction to constant w allows one to discard considerations related 
to specific interferometer configurations and to study the spread spectrum 
phenomenon at light computational load \6j. 



where the vector n represents noise corrupting the measure- 
ments. The M X N measurement matrix $ identifies the 
complete linear relation between the signal and the visibili- 
ties. The matrix A is the diagonal matrix implementing the 
primary beam. The matrix C is the diagonal matrix imple- 
menting the modulation by the -modulation, which reduces 
to the identity if small FoV assumptions are made. The uni- 
tary matrix F implements the discrete Fourier transform. The 
M xN matrix M is the rectangular binary matrix implement- 
ing the mask that characterises visibility coverage. Given in- 
complete visibility coverage with M < N, ^ defines an 
ill-posed inverse problem. 

3. COMPRESSED SENSING 
FOR INTERFEROMETRIC IMAGING 

CS reconstruction techniques have been applied to interfer- 
ometric imaging recently |5|, outperforming the standard 
CLEAN algorithm. The spread spectrum phenomenon which 
arises naturally in radio interferometry, provides a further en- 
hancement in reconstruction quality |6|. Extensions to wide 
field imaging provide another enhancement |7|. We review 
these works in this section. 

3.1. Sparsity for interferometric image reconstruction 

The theory of CS acknowledges that a large variety of sig- 
nals in Nature are sparse. In this context, the interferometric 
ill-posed inverse problem can be regularised by a sparsity or 
compressibility prior. This results in an optimisation problem 
called Basis Pursuit (BP), which consists of minimising the 
^i-norm of the coefficients a of the signal represented in a 
sparsity basis ^ (the N x N matrix with columns containing 
the vectors of the sparsity basis), under a constraint on the 
^2 -norm of the residual noise: 

rnin||a||i such that ||y - <l>^a||2 < e . 

Recall that the £i-norm is given by ||a||i = Y^f^i while 
the -^2-norm is given by ||a||2 = (Xl^^i |<^iP)"^^^. The con- 
straint e may be related to a residual noise level estimator fSl. 
The BP problem may be solved by the application of non- 
linear, iterative convex optimisation algorithms. If the solu- 
tion of the optimisation problem is denoted a^, then the sig- 
nal is reconstructed through the synthesis x'^ = ^a^. 

Under typical small FoV assumptions, solving the BP 
problem provides similar reconstruction quality to multi- scale 
versions of the CLEAN algorithm. However, the versatility 
of this framework allows the easy addition of many useful 
priors, such as the positivity prior leading to BP+ reconstruc- 
tions of superior quality to CLEAN |5 1. A typical simulation 
used to verify this result and the reconstruction quality of var- 
ious imaging techniques, as a function of the percentage of 
visibilities measured, are shown in Fig. [2] The Dirac sparsity 
basis was used in this work. 




Fig. 2. BP, BP+ and CLEAN reconstruction performance on 
elongated Gaussian profile signals (credit: f5\). 

3.2. Spread spectrum phenomenon 

If small FoV assumptions are not made, then the -modulation 
modulation C no longer reduces to the identity. Beyond 
sparsity, CS demonstrates that the incoherence between the 
measurement and sparsity bases is a key criterion driving the 
performance of CS based reconstructions. In the interfero- 
metric setting the measurement basis is essentially the Fourier 
basis, hence coherence is given by the maximum modulus of 
the Fourier coefficient of the sparsity basis vectors. Con- 
sequently, an operation that acts to reduce the maximum 
Fourier coefficient, reduces the coherence and thus improves 
the quality of CS reconstructions. The it; -modulation cor- 
responds to a norm-preserving convolution in the Fourier 
plane, spreading the spectrum of the sparsity basis vectors, 
and achieves exactly this. The Dirac basis considered pre- 
viously is already maximally incoherent, hence this spread 
spectrum phenomenon is not expected to improve perfor- 
mance for this case. However, for alternative bases we expect 
an enhancement. 

Many natural signals are also sparse or compressible in 
the magnitude of their gradient, in which case the Total Vari- 
ation (TV) minimisation problem applies |3|. The TV prob- 
lem involves replacing the £i-norm sparsity prior in the BP 
problem with a prior on the TV norm of the signal itself: 

rmn||a?||TV such that \\y — ^x\\2 < e ^ 

where the TV norm is defined by the ^i-norm of the gradi- 
ent of the signal ||^c||tv = ||Vx||i. The TV problem may 
also be solved by the application of non-linear, iterative con- 
vex optimisation algorithms, illustrating another example of 
the versatility of the framework. The signal is directly re- 
covered from the solution to the optimisation problem. Re- 
covering interferometric images by solving the TV problem 
has been shown to achieve very good reconstruction quality 
and, moreover, is enhanced considerably by the spread spec- 
trum phenomenon |6|. A typical simulation used to study this 
effect and reconstruction quality for BP (with Dirac sparsity) 



Fig. 3. BP and TV reconstruction performance on Gaussian 
profile signals, in the absence and presence of the spread spec- 
trum (SS) phenomenon (credit: L6J). 

and TV approaches, in the absence and presence of the spread 
spectrum phenomenon, are shown in Fig. [3] 

3.3. Wide field-of-view 

Incorporating wide FoV considerations in interferometric 
imaging is of increasing importance for the next generation 
of radio interferometric telescopes. CS techniques for in- 
terferometric imaging have been extended to a wide FoV 
recently |7|, recovering images in the spherical coordinate 
space in which they naturally live. The effectiveness of the 
spread spectrum phenomenon is enhanced when going to a 
wide FoV since the it; -modulation evolves to include higher 
frequency content, while sparsity is promoted by recovering 
images directly on the sphere. Both of these properties act to 
improve the quality of reconstructed interferometric images. 
The wide FoV ill-posed inverse problem is defined by 

y ^sXs + n , (3) 

with 

= MFCAP, 

where the vector Xs is the interferometric image defined on 
the celestial sphere (rather than a tangent plane). The mea- 
surement operator on the sphere simply consists of aug- 
menting the operator on the plane $ by prepending a projec- 
tion P from the sphere to the plane, which incorporates a con- 
volutional regridding. Careful consideration is given to sam- 
plings on the sphere and plane to ensure that the planar grid is 
sampled sufficiently to accurately represent the projection of 
a band-limited signal defined on the sphere. The BP and TV 
problems are then solved to recover interferometric images 
on the sphere. This approach has been shown to perform very 
well, with the quality of recovered spherical interferometric 
images superior to planar reconstructions |7|. A typical sim- 
ulation used to study this effect and reconstruction quality for 
BP (with Dirac sparsity) and TV approaches both on the pla- 
nar and spherical settings, in the absence and presence of the 
spread spectrum phenomenon, are shown in Fig.|4] 
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(a) Typical simulation 



(b) BP reconstruction performance 



(c) TV reconstruction performance 



Fig. 4. BP and TV reconstruction performance on Gaussian profile signals on the plane and sphere, in the absence and presence 
of the spread spectrum (SS) phenomenon (credit: (author?) | 7 1). 



4. FUTURE DIRECTION 

The CS techniques developed for interferometric imaging 
E m [3 remain somewhat idealised. Random visibility 
coverage in (li, v) is assumed, with the spatial frequencies 
probed by the interferometer also assumed to fall on discrete 
grid points. Furthermore, to study the spread spectrum phe- 
nomenon a constant w is assumed. These restrictions have 
been necessary to remain as close to the theory of CS as possi- 
ble during the development and evaluation of interferometric 
imaging techniques. However, now that the effectiveness of 
these techniques has been demonstrated, it is important to 
adapt them to realistic interferometric configurations. 

In practice a varying w will not alter the essence of the 
spread spectrum phenomenon due to its realisation through 
the convolution operator; however, this effect must be studied. 
The i(;-projection algorithm | 9 1 may be applied to reduce the 
computational cost in this context. 

Realistic visibility coverages in all components {u, w) 
must also be considered by simulating observations made by 
real interferometer configurations. Furthermore, to integrate 
CS imaging techniques into the data processing pipelines of 
radio interferometric telescopes, it is necessary to consider the 
continuous visibilities measured by an interferometer before 
they are gridded (so that imaging can be integrated into the 
iterative self-calibration of telescopes). 

In general, CS addresses imaging by optimising both re- 
construction and acquisition, while current applications in in- 
terferometry have essentially focused on reconstruction only. 
The possibility of optimising the configuration of interferom- 
eters to enhance the spread spectrum phenomenon 1 6 1 for CS 
reconstruction is an exciting avenue of research at the level of 
acquisition, contrary to the traditional approach that attempts 
to minimise this effect. In addition, direction dependent beam 
effects may also provide an alternative source of the spread 
spectrum phenomenon. 

In conclusion, with a new generation of radio interfero- 
metric telescopes under design and construction, such as the 



Square Kilometer Array (SKA), the development of novel in- 
terferometric techniques is of paramount importance. It is es- 
sential that the CS interferometric imaging techniques devel- 
oped recently are adapted to realistic settings to ensure that 
the fidelity of reconstructed images keeps pace with the capa- 
bilities of new instruments. 
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